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Abstract 


Prompt photon production in ^JsNN - 2.76 TeV Pb+Pb collisions has been measured by the 
ATLAS experiment at the Large Hadron Collider using data collected in 2011 with an integ¬ 
rated luminosity of 0.14 nb“^ Inclusive photon yields, scaled by the mean nuclear thickness 
function, are presented as a function of collision centrality and transverse momentum in two 
pseudorapidity intervals, |?7| < 1.37 and 1.52 < \r]\ < 2.37. The scaled yields in the two pseu¬ 
dorapidity intervals, as well as the ratios of the forward yields to those at midrapidity, are 
compared to the expectations from next-to-leading order perturbative QCD calculations. The 
measured cross sections agree well with the predictions for proton-proton collisions within 
statistical and systematic uncertainties. Both the yields and ratios are also compared to two 
other pQCD calculations, one which uses the isospin content appropriate to colliding lead 
nuclei, and another which includes nuclear modifications to the nucleon parton distribution 
functions. 
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1 Introduction 

Prompt photons are an important probe for the study of the hot, dense matter formed in the high-energy 
collision of heavy ions. Being colorless, they are transparent to the subsequent evolution of the matter 
and probe the very initial stages of the collision. Their production rates are therefore expected to be 
directly sensitive to the overall thickness of the colliding nuclear matter. The rates are also expected to be 
sensitive to modifications of the partonic structure of nucleons bound in a nucleus, which are implemented 
as nuclear modifications [1, 2, 3] to the parton distribution functions (PDF) measured in deep-inelastic 
lepton-proton and proton-proton (pp) scattering experiments. These effects include nuclear shadowing 
(the depletion of the parton densities at low Bjorken x), anti-shadowing (an enhancement at moderate x), 
and the EMC effect [4]. Photon rates are also sensitive to final-state interactions in the hot and dense 
medium, via the conversion of high energy quarks and gluons into photons through rescattering. This is 
predicted to lead to an increased photon production rate relative to standard expectations [5, 6]. 

Prompt photons have two primary sources. The first is direct emission, which proceeds at leading order 
via quark-gluon Compton scattering qg qy or quark-antiquark annihilation qq gy. The second is the 
fragmentation contribution from the production of hard photons during parton fragmentation. At leading 
order in perturbative quantum chromodynamics (pQCD) calculations, there is a meaningful distinction 
between the direct emission and fragmentation, but at higher orders the two cannot be unambiguously 
separated. In order to suppress the large background of nonprompt photons originating from the decays 
of neutral mesons in jets, as well as fragmentation photons, an isolation criterion is applied, both in meas¬ 
urements and calculations, to the transverse energy contained within a cone of well-defined size around 
the photon direction [7]. The isolation transverse energy requirement can be applied as a fraction of the 
photon transverse energy, or as a constant transverse energy threshold. In either case, these requirements 
can be applied consistently to pQCD calculations so that prompt photon rates can be calculated reliably, as 
the isolation criterion naturally cuts off the collinear divergence of the fragmentation contribution [7]. 
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Prompt photon rates have been measured extensively in both fixed target and collider experiments. Fixed 
target experiments include WA70 [8], UA6 [9] and E706 [10], and cover the range ^/s = 23-38.8 GeV. In 
collider experiments, measurements were performed for proton-proton collisions at the CERN Intersect¬ 
ing Storage Rings {pp, yfs = 24-62.4 GeV ) [11, 12], and BNE Relativistic Heavy Ion Collider {pp at 
yfs =200 GeV) [13, 14], and for proton-antiproton collisions at the CERN Super Proton Synchrotron {pp, 
yfs = 546-630 GeV) [15, 16] and at the Eermilab Tevatron {pp, - 0.63-1.96 TeV) [17, 18, 19, 20]. 
At the CERN Earge Hadron Collider (EHC), ATEAS [21, 22, 23] and CMS [24, 25] have measured 
isolated prompt photons in pp collisions at = 7 TeV. In most cases, good agreement has been found 
with pQCD predictions at next-to-leading order (NEC), which are typically calculated using the Jetphox 
package [7, 26]. In lower-energy heavy-ion collisions, the WA98 experiment observed direct photons in 
lead-lead (Pb-i-Pb) collisions at a/^nn - 17.3 GeV [27], and the PHENIX experiment performed meas¬ 
urements of direct photon rates in gold-gold collisions at a/Snn - 200 GeV [28, 29]. 

A variable often used to characterize the modification of rates of hard processes in a nuclear environment 
is the nuclear modification factor 

_ {llN,,ddNxldpT 

{TAA}d(rPx"ldpT ’ ^ 

where dNxIdpT is the yield of objects X produced in a px interval, Vgvt is the number of sampled 
minimum-bias events, Taa is the mean nuclear thickness function (defined as the mean number of binary 
collisions divided by the total inelastic nucleon-nucleon {NN) cross section) and dcr'^/dpj is the cross 
section of process X in pp collisions for the same px interval. With this formula, one can make straight¬ 
forward comparisons of yields in heavy-ion collisions, normalized by the flux of initial-state partons, to 
those measured in pp data, or calculated in pQCD. CMS performed the first measurement of isolated 
prompt photon rates in both Pb-i-Pb and pp collisions at ^/s =2.76 TeV up to a photon transverse energy 
Ex = 80 GeV within \r]\ < 1.44 [30]. This measurement observed prompt, isolated photon rates consistent 
with Raa - 1 for all collision impact parameters and Ej ranges considered, and good agreement of the 
data with Jetphox calculations. 

This paper presents isolated prompt photon yields, scaled by the mean nuclear thickness to derive ef¬ 
fective cross sections, measured in Pb-i-Pb collisions with the ATEAS detector, making use of its large- 
acceptance, longitudinally segmented calorimeter system. The effect of the underlying event (UE) on the 
photon energy and shower shape is corrected on an event-by-event basis. Photon yields are measured 
over two ranges in the pseudorapidity of the photon, Ipl < 1.37 (central) and 1.52 < Ipl < 2.37 (forward), 
and for photon transverse momenta in the interval 22 < px < 280 GeV. Comparisons of the yields with 
NEO pQCD calculations are also presented from Jetphox 1.3 [26], in three configurations: pp collisions, 
Pb-i-Pb collisions (i.e. with the correct total isospin), and Pb-i-Pb after incorporating the EPS09 nuclear 
modification factors to the nucleon PDEs [1], derived from experimental data of lepton and proton scat¬ 
tering on nuclei. The ratios of the yields in the forward rj region to those in the central rj region (RpC;?) are 
also presented. 


2 ATLAS detector 

The ATEAS detector comprises three major subsystems: the inner detector, the calorimeter system, and 
the muon spectrometer. It is described in detail in Ref. [31]. 
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The inner detector is comprised of the pixel detector, the semiconductor tracker (SCT) and the transition 
radiation tracker (TRT), which cover the full azimuthal range and pseudorapidities ^ |?7| < 2.5, except for 
the TRT, which covers \t]\ < 2. The muon spectrometer measures muons over I77I < 2.7 with a combination 
of monitored drift tubes and cathode strip chambers. 

The ATLAS calorimeter is the primary subsystem used for the measurement presented here. It is a large- 
acceptance, longitudinally segmented sampling calorimeter covering I77I < 4.9 with electromagnetic (EM) 
and hadronic sections. The EM section is a lead/liquid-argon sampling calorimeter with an accordion¬ 
shaped geometry. It is divided into a barrel region, covering \r]\ < 1.475, and two endcap regions, covering 
1.375 < I77I < 3.2. The EM calorimeter has three primary sections, longitudinal in shower depth, called 
“layers,” to fully contain photon showers in the range of interest for this analysis. The first sampling 
layer is 3 to 5 radiation lengths deep and is segmented into fine strips of size Ary = 0.003 - 0.006 
(depending on rf), which allows the discrimination of photons from the two-photon decays of 7r° and 
T] mesons. The second layer is 17 radiation lengths thick, sampling most of an electromagnetic shower, 
and has cells of size Arjx A(p = 0.025 x 0.025. The third layer has a material depth ranging from 4 to 
15 radiation lengths and is used to correct for the leakage beyond the first two layers for high-energy 
electromagnetic showers. The total material in front of the electromagnetic calorimeter ranges from 2.5 
to 6 radiation lengths depending on pseudorapidity, except in the transition region between the barrel and 
endcap regions (1.37 < |?7| < 1.52), in which the material is up to 11.5 radiation lengths (for which reason 
this transition region is excluded from this analysis). In front of the strip layer, a presampler is used to 
correct for energy loss in front of the calorimeter within the region |?7| < 1.8. In test beam environments 
and in typical pp collisions, the photon energy resolution is found to have a sampling term of 10-17% 
/ V-EfGeV]. Above 200 GeV the global constant term in the photon energy resolution, estimated to be 
1.2%+0.6% (1.8%±0.6%) in the barrel (endcap) region for ppdataat ^/s = 7TeV, starts to dominate [32]. 
The hadronic calorimeter section is located outside the electromagnetic calorimeter. Within \r]\ < 1.7, it 
is a sampling calorimeter of steel and scintillator tiles, with a depth of 7.4 hadronic interaction lengths. 

The ATEAS zero-degree calorimeters (ZDCs) are used for minimum-bias event triggering. They detect 
forward-going neutral particles with I77I > 8.3. The minimum-bias trigger scintillators (MBTS) detect 
charged particles in the interval 2.1 < \t]\ < 3.9 using two sets of 16 counters positioned at z = ±3.6 m. 
They are used for event selection. The forward calorimeter (ECal) is used to determine the “centrality” of 
the collision, which can be related to geometric parameters such as the number of participating nucleons 
or the number of binary collisions [33]. The ECal has three layers in the longitudinal direction, one 
electromagnetic and two hadronic, covering 3.1 < |?y| < 4.9. The ECal electromagnetic and hadronic 
modules are composed of copper and tungsten absorbers, respectively, with liquid argon as the active 
medium, which together provide ten interaction lengths of material. 

The sample of events used in this analysis was collected using the first-level calorimeter trigger [34]. 
This is a hardware trigger that sums the electromagnetic energy in towers of size A/y x A0 = 0.1 x 0.1. A 
sliding window of size 0.2 x 0.2 was used to find elecfromagnefic clusfers by searching for local energy 
maxima and keeping only fhose clusters wifh energy in fwo adjacenf lowers (i.e. regions wifh a size of 
eilher 0.2 x 0.1 or 0.1 x 0.2) exceeding a Ihreshold. The trigger used for the present measurement had a 
threshold of 16 GeV transverse energy. 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the y axis points 
upward. Cylindrical coordinates (r, (p) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(0/2). 
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3 Collision data selection 


The data sample analyzed in this paper corresponds to an integrated luminosity of Xint = 0.14 nb“^ 
Pb+Pb collisions at ^JsNN = 2.76 TeV collected during the 2011 LHC heavy-ion run. After the trigger 
requirement, events must satisfy a set of selection criteria. To suppress backgrounds, the relative time 
measured between the two MBTS counters is required to be less than 5 ns, and a primary vertex is 
required to be reconstructed in the inner detector. Minimum-bias events were triggered in the same data 
samples based on either a coincidence in the two ZDCs associated with a track in the inner detector, or a 
total of at least 50 GeV transverse energy deposited in the full calorimeter system. These events were also 
required to pass the same MBTS and vertex selections as the photon-triggered events. To be consistent 
with the minimum-bias trigger selections, a ZDC coincidence is also required for photon-triggered events 
with low FCal XXt- 

The centrality of each heavy-ion collision is determined using the total transverse energy measured in 
the forward calorimeter (3.2 < I 77 I < 4.9), at the electromagnetic scale, FCal 'LEj. The trigger and event 
selection were studied in detail in the 2010 Pb-i-Pb data sample [35] and 98 ±2% of the total inelastic cross 
section was accepted. The higher luminosity of the 2011 heavy-ion run necessitated a more sophisticated 
trigger strategy, including more restrictive triggers in the most peripheral events. However, it was found 
that the FCal EXj distributions in 2011 data match those measured in 2010 to a high degree of precision. 
For this analysis, the FCal l,Ej distribution was divided into four centrality intervals, covering the 0- 
10%, 10-20%, 20^0% and 40-80% most central events. With this convention, the 0-10% interval 
contains the events with the largest forward transverse energy production, and the 40-80% interval the 
smallest. The total number of minimum-bias events corresponding to the 0-80% centrality interval is 
Aevt - 7.93 X 10^. 

Quantities which describe the average geometric configuration of the colliding nuclei are calculated as 
described in Ref. [36] using a Glauber Monte-Carlo calculation to describe the measured minimum-bias 
FCal distribution. Table 1 summarizes all of the centrality-related information used in this analysis. For 
each centrality interval, the table specifies fhe mean number of nucleons fhaf inferacf af leasf once (A^part)^ 
fhe mean number of binary collisions (77coii), and fhe mean value of fhe nuclear fhickness function (Taa), 
wifh fheir respecfive fraclional uncerfainfies. The uncertainly on fhe mean nuclear fhickness function 
(Taa) = {Ncon)/o'NN is smaller fhan fhe corresponding uncertainly on (77coll)^ since fhe uncertainly on 
(Tnn largely cancels in fhe ralio. All of fhe uncerfainfies accounl for variafions in fhe Glauber model 
parameters consisfenf wifh fhe uncerfainfies aboul fhe nuclear wave funclion, as well as fhe uncerfainfy in 
fhe estimation of fhe measured fracfion of fhe lolal inelaslic cross seclion. 

Since fhe dislribulion of FCal EFt is differenl in evenls wifh high-pj phofons compared lo minimum-bias 
evenls, a weighfing faclor is applied fo each simulafed even! lo make fhe simulaled dislribufions agree 
wifh fhe measured dislribufions. 

4 Simulated data samples 

For fhe exlraclion of photon reconsfruclion and idenlificalion efficiencies, fhe pholon energy scale, and 
expecled properties of fhe isolation Iransverse energy dislribufions, samples of evenls confaining prompl 
photons were produced using Pythia 6.423 [37] for pp collisions at ^/s = 2.76 TeV using fhe ATLAS 
AUET2B sel of luned paramelers [38]. Direcl phofons were simulafed in pholon-jel evenls divided info 
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Table 1: Centrality bins used in this analysis, tabulating the percentage range, the average number of participants 
{{Npsn}) and binary collisions ((A^coii)), the mean nuclear thickness ((Taa)) and the relative systematic uncertainty 
on these quantities. 


Interval 

(Apart) 

^(^part) 

(Apart) 

(Acoll) 

5<Aco11> 

(Acoll) 

(Taa) [mb-i] 

6{Tkk) 

(Taa) 

0-10% 

356.2 

0.7% 

1500.6 

7.6% 

23.4 

1.6% 

10-20% 

260.7 

1.4% 

923.3 

7.3% 

14.4 

2.1% 

20-40% 

157.8 

2.4% 

440.6 

7.2% 

6.9 

3.5% 

40-80% 

45.9 

5.9% 

77.8 

9.1% 

1.2 

8.1% 


four sub-samples based on requiring a minimum pj for the primary photon: pj > 17 GeV, pj > 35 GeV, 
Pt > 70 GeV, and pT > 140 GeV. The contribution of fragmentation photons was modeled using a set 
of simulated inclusive jet pp events, also using the same Pythia 6 tune. Each of these is required to 
have a hard photon produced in the fragmentation of jets produced with the Pythia 6 hardness scale, 
which controls the typical px of the produced jets, ranging from 17 to 560 GeV. Similar samples were 
also prepared using the Sherpa generator [39] using the CTIO [40] parton distribution functions, which 
include both direct and fragmentation photon contributions. These were used to check on the generator 
dependence of the photon efficiency. A large sample of Pythia 6 inclusive-jet events, without the hard 
photon requirement, were utilized to study the properties of background candidates. For all generated 
samples, each event was fully simulated using GEANT4 [41, 42]. 

Each simulated event is overlaid upon a real minimum-bias event from experimental data, with the simu¬ 
lated event vertex placed at the position of the measured vertex position. By using minimum-bias data as 
the underlying-event model, almost all features of the underlying event are preserved in the simulation, 
including the full details of its azimuthal correlations. 

A reconstructed photon is considered “matched” to a prompt generator-level (“truth”) photon when they 
are separated by an angular distance AR = < 0.2. If multiple reconstructed photons are 

within the matching window, only the highest-px reconstructed photon is considered matched to the truth 
photon. 


5 Photon reconstruction 

The electromagnetic shower associated with each photon, as well as the total transverse energy in a 
cone surrounding it, are reconstructed as described in Ref. [43]. However, in a heavy-ion collision, it is 
important to subtract the large UE from each event before the reconstruction procedure is applied. If it is 
not subtracted, photon transverse energies can be overestimated by up to several GeV in the most central 
events, and the isolation transverse energy in a AR = 0.3 cone can be overestimated by about 60 GeV. 
The procedure explained in Ref. [44] is used to estimate the energy density of the underlying event in 
each calorimeter cell. It iteratively excludes jets from consideration in order to obtain the average energy 
density in each calorimeter layer in intervals of Arj = 0.1, after accounting for the elliptic modulation 
relative to the event plane angle measured in the FCal [35, 45]. The algorithm provides the energy density 
as a function of p, (p, and calorimeter layer, which allows the event-by-event subtraction of the UE in the 
electromagnetic and hadronic calorimeters. 
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After subtraction, the residual deposited energies stem primarily from three sources: jets, photons/electrons, 
and UE fluctuations (including higher-order flow harmonics). It should be noted that while this provides 
an estimate of the mean underlying transverse energy as a function of 77 , it is at present not possible to 
make further subtraction of more localized structures. 

The ATLAS photon reconstruction [43] is seeded by clusters with £7 > 2.5 GeV found using a sliding- 
window algorithm applied to the second sampling layer of the electromagnetic calorimeter, which typ¬ 
ically contains over 50% of the shower energy. In the dense environment of the heavy-ion collision, the 
photon conversion reconstruction procedure is not performed, due to the large number of combinatoric 
pairs in more central collisions. However, a substantial fraction of converted photons are still reconstruc¬ 
ted by the photon algorithm as, for high energy photon conversions, the electron and positron are typically 
close together when they reach the calorimeter, while their tracks typically originate at a radius too large 
to be well described by the tracking algorithm that is used for heavy-ion collisions. Thus, the photon 
sample analyzed here is a mix of converted and unconverted photons. From simulations, the overall 
conversion rate is found to be about 30% in |? 7 | < 1.37 and 60% in 1.52 < |? 7 | < 2.37. 

The energy measurement is made using the three layers of the electromagnetic calorimeter and the pres¬ 
ampler, with a window size corresponding to 3 x 5 cells (in rj and cfi) in the second layer in the barrel, 
and 5x5 cells in the endcap region. An energy calibration is applied to each shower to account for both 
its lateral leakage (outside the nominal window) and longitudinal leakage (into the hadronic calorimeter 
as well as dead material) [43]. For converted photons, this window size can lead to an underestimate of 
the photon candidate’s energy, which is accounted for in the data analysis. The transverse energy of the 
photon is defined as the calibrated cluster energy multiplied by the sine of the polar angle determined 
with respect to the measured event vertex. The transverse momentum of the photon is identified with the 
measured transverse energy. 

The fine-grained, longitudinally segmented calorimeter allows for a detailed characterization of the shape 
of each photon shower, which can be used to reject neutral hadrons while maintaining a high efficiency 
for photons. Nine shower shape variables are used for each photon candidate: 

The primary shape variables used can be broadly classified by which sampling layer is used. The second 
sampling layer is used to measure 

• the ratio of energies deposited in a 3 x 7 (77 x 0) window to those deposited in a 7 x 7 set of 
cells in the second layer. 

• Rf. the ratio of energies deposited in a 3 x 3 (77 x 0) window to those deposited in a 3 x 7 set of 
cells in the second layer. 

• WrjX- the standard deviation in the 77 projection of the energy distribution of the cluster in a 3 x 5 set 
of cells in the second layer. 

The hadronic calorimeter is used to measure the fraction of shower energy that is detected behind the 
electromagnetic calorimeter. Only one of these is applied to each photon, depending on its pseudorapid¬ 
ity: 

• ^had^ the ratio of transverse energy measured in the hadronic calorimeter to the transverse energy 
of the photon candidate. This quantity is used for 0.8 < I 77 I < 1.37. 

• ^hadi: the ratio of transverse energy measured in the first sampling layer of the hadronic calorimeter 
to the transverse energy of the photon candidate. This quantity is used for photons with either 
I 77 I < 0.8 or I 77 I > 1.52. 
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Finally, cuts are applied to five other quantities, measured in the fine-granularity first layer, to reject 
neutral meson decays from jets. In this finely segmented layer a search for multiple maxima from elec¬ 
tromagnetic decays of neutral hadrons is performed: 

• Ws^tot- the standard deviation of the energy distribution in the 77 projection in the first sampling 
“strip” layer, in strip cell units. 

• : the standard deviation of the energy distribution in three strips including and surrounding the 
cluster maximum in the strip layer, also in strip cell units. 

• ^side^ the fraction of energy in seven strips surrounding the cluster maximum, not contained in the 
three core strips. 

• ^^ratio^ the asymmetry between the energies in the first and second maxima in the strip layer cells. 
This quantity is equal to one when there is no second maximum. 

• AE: the difference between the energy of the second maximum, and the minimum cell energy 
between the first two maxima. This quantity is equal to zero when there is no second maximum. 

In a previous ATLAS measurement [21], it was observed that the distributions of the shower shape vari¬ 
ables measured in data differ systematically from those in the simulation. To account for these differ¬ 
ences, a set of correction factors was derived, each of which changes the value of a simulated shower 
shape variable such that its mean value matches that of the corresponding measured distribution. For the 
measurements presented in this paper, the same correction factors, obtained by comparing pp simulations 
to the same quantities in data, are used with no modification for the heavy-ion environment. They were 
validated in the heavy-ion environment using electrons and positrons from reconstructed Z ^ e'^e~ de¬ 
cays from the same LHC run. It was observed that the magnitude and centrality dependence of the mean 
values of the shape variables are well described by simulations, within the limited size of the electron and 
positron sample. 

Figure 1 shows three typical distributions of shower shape variables for data from the 0-10% and 40-80% 
centrality intervals, each compared with the corresponding quantities in the simulation. The simulated 
distributions, after shower shape corrections, are all normalized to the number of counts in the corres¬ 
ponding data histogram. The data contain some admixture of neutral hadrons, so complete agreement 
should not be expected in the full distributions. The admixture of converted photons, which depends on 
the amount of material in front of the electromagnetic calorimeter, and thus the pseudorapidity of the 
photon, is not accounted for in the analysis, but there is good agreement of the shower shape variable dis¬ 
tributions between data and simulation. Converted photons tend to have wider showers than unconverted 
photons, and so substantially broaden the shower shape variables. 

The electromagnetic-energy-trigger efficiency was investigated using a sample of minimum-bias data, 
where the primary triggers did not select on particular high-px activity. Using these, the probability 
for photon candidates selected for this analysis to match a first-level trigger with F'x.trig >16 GeV and 
AR < 0.15 exceeds 99% for well-reconstructed photon candidates with px ^ 22 GeV and over the full 
centrality range. In the more central events, the underlying-event contribution to the photon candidate 
reduces the effective threshold down by several GeV relative to the more peripheral events. To work in 
the plateau region, the minimum px required in this analysis is 22 GeV. 

Photons are selected for offline analysis using a variation of the “tight” selection criteria developed for 
the photon analysis in pp collisions [ 21 ], necessitated by the additional fluctuations in the shower shape 
variables induced by the underlying event in heavy-ion collisions. Specific intervals are defined for all 
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Figure 1; (Color online) Comparisons of distributions of three shower shape variables (ws, 3 , and f?had) from 
data (black points) with simulation results after shower shape corrections (yellow histogram), for tight and isolated 
photons with reconstructed 35 < px < 44.1 GeV and I 77 I < 1.37. Events from the 0-10% centrality interval are 
shown in the top row (a)-(c), while those from the 40-80% interval are shown in the bottom row (d)-(f). 


nine shower shape variables, and are implemented in a px-independent, but //-dependent scheme. The 
intervals for each variable are defined to contain 97% of the distribution of isolated reconstructed photons 
matched to isolated truth photons with a reconstructed px in the region 40 < px < 60 GeV in the 0-10% 
centrality interval (where the UE fluctuations are largest), using the isolation criteria described in the next 
section. 

In order to derive a data-driven estimate of the background candidates from jets, a “nontight” selection 
criterion is defined, which is particularly sensitive to neutral hadron decays. For this selection, a photon 
candidate is required to fail at least one of four shower shape selections in the first calorimeter layer: 

Eside> Eratio and AE. These reversed selections enhance the probability of accepting neutral hadron decays 
from jets, via candidates with a clear double shower structure (via Eratio and AE) as well as candidates in 
which the two showers may have merged (via 11 ) 5,3 and Eside) [ 21 ]. 

While the photon energy calibration is the same as used for pp collisions, based in part on measurements 
of Z bosons decaying into an electron and a positron, and validated with Z ^ {i + y events [46], the 
admixture of converted and unconverted photons leads on average to a small underestimate of the photon 
energy in Pb-i-Pb events, since the energies of converted photon clusters is typically reconstructed in a 
larger region in the calorimeter. This is quantified in the simulation by the mean fractional difference 
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between the reconstructed and the truth photon transverse momenta 

obtained from simulation. For matched photons, the average deviation from the truth photon pj is the 
largest at low photon pj and is typically within 1 % for px > 44 GeV. The fractional energy resolution, 
determined by calculating the standard deviation of the same quantity in smaller intervals in ranges 
from 4.5% for 22 < p™*’ < 26 GeV to 1.5% for p™*’ = 200 GeV for \rj\ < 1.37 and from 6% to 3% 
for 1.52 < \t]\ < 2.37. The effects of energy scale and resolution are corrected for by using bin-by-bin 
correction factors described below. 

The isolation transverse energy is the sum of transverse energies in calorimeter cells (including 
hadronic and electromagnetic sections) in a cone of size ARiso around the photon direction. The photon 
energy is removed by excluding a central core of cells in a region corresponding to 5x7 cells in the second 
layer of the EM calorimeter. The cone size is chosen to be A7?iso - 0.3, to reduce the sensitivity to UE 
fluctuations. The isolation criterion is E™ < 6 GeV. An additional correction, based on simulations and 
parameterized primarily by the photon energy and p, is then applied to the calculated isolation transverse 
energy to minimize the effects of photon shower leakage into the isolation cone. It typically amounts to a 
few percent of the reconstructed photon transverse energy. 

The left column of Eig. 2 shows the distributions of E™ for tight photon candidates with 35 < px < 
44.1 GeV as a function of collision centrality, compared with simulated distributions. The data and 
simulations are normalized so the integrals of E}^° < 0, where no significant background from jet events 
is expected, are the same. Both, the simulated and the measured Ej° distributions grow noticeably wider 
with increasing centrality; as the UE subtraction only accounts for the mean energy in an rj interval, 
local fluctuations are still present. Eurthermore, in the data, an enhancement in events with E™ > 0 is 
expected from the jet background. The Ej° distribution for a sample enhanced in backgrounds is shown 
in the right column of Eig. 2, which shows the isolation distribution for the nontight candidates in the 
same px interval. Eor larger values of Ej°, the distributions from the tight and nontight samples have 
similar shapes. The distributions are normalized to the integral of the tight photon candidate distribution 
in the region Ej° > 8 GeV. 

After applying the tight selection and an isolation criterion of E}^° < 6 GeV to the 0-80% centrality 
sample, there are 62 130 candidates with px > 22.0 GeV within |p| < 1.37 and 30568 candidates within 
1.52 < IpI < 2.37. 


6 Yield extraction 

The kinematic intervals used in this analysis are defined as follows. Eor each centrality interval, as 
described in Sect. 3, the photon kinematic phase space is divided into intervals in photon p and pj. The 
two primary regions in p are |p| < 1.37 (“central p”), 1-52 < |p| < 2.37 (“forward p”). The px intervals 
used are logarithmic, and are 17.5 < px < 22 GeV (only used in simulations), 22.0 < px < 27.8 GeV, 
27.8 < Px < 35.0 GeV, 35.0 < px < 44.1 GeV, 44.1 < px < 55.6 GeV, 55.6 < px < 70.0 GeV, 
70.0 < Px < 88.2 GeV, 88.2 < px < 140 GeV, and 140 < px < 280 GeV. 

Prompt photons are defined as photons produced in the simulation of the hard process, either directly or 
radiated from a primary parton, via a truth particle-level isolation transverse energy selection of Ej° < 
6 GeV. The truth-level Ej'^ is defined using all final-state particles except for muons and neutrinos in a 
cone of AR = 0.3 around the photon direction. To account for the underlying event in the hard process. 
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Figure 2; (Color online) Distributions of photon isolation transverse energy in a Af^iso = 0.3 cone for the four 
centrality bins in data (black points, normalized by the number of events and by the histogram bin width), for 
photons with 35 < px < 44.1 GeV. In the left column (a)-(d) simulations (yellow histogram) are normalized to 
the data so that the integrals in the range Ej° < 0 are the same. The corresponding sample of nontight photon 
candidates, normalized to the distribution of tight photons for is™ > 8 GeV is shown overlaid on the tight photon 
data in the right column (e)-(h) to illustrate the source of the photons with large Ej°. 


the mean energy density is estimated for each simulated event using the jet-area method described in 
Ref. [21]. 

For each interval in pj, rj and centrality (C), the per-event yield of photons is defined as 

1 dNy Nfai{pT,p,cw{pT,n,c) 

fVevt(C)t/pT ” (Vevt(C)6tot(PT,t7,C)ApT ’ 

where is the background-subtracted yield, 'W is a factor that corrects for the bin migration due to the 
photon energy resolution and any residual bias in the photon-energy scale, 'TF is a factor that corrects for 
electron contamination from W and Z bosons, etot is the combined photon reconstruction and identification 
efficiency, A^evt is the number of minimum-bias events in centrality interval C, and Apx is the width of the 
transverse momentum interval. 

The technique used to subtract the background from jets from the measured yield of photon candidates 
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Figure 3: (Color online) Illustration of the double sideband approach, showing the two axes for partitioning photon 
candidates: region A is the “signal region” (tight and isolated photons), region B contains tight, nonisolated photons, 
region C contains nontight isolated photons, and region D contains nontight and nonisolated photons. 


is the “double sideband” method, used in Refs. [21, 22, 23]. In this method, photon candidates are 
partitioned along two dimensions, illustrated in Fig. 3. The four regions are labeled A, B, C and D and 
correspond to the four categories expected for reconstructed photons and background candidates: 

• A: tight, isolated photons: Signal region for prompt, isolated photons. 

• B: tight, nonisolated photons: A region expected to contain nonisolated photons produced in the 
vicinity of a jet or an upward UE fluctuation, as well as hadrons from jets with shower shapes 
similar to those of a tight photon. 

• C: nontight, isolated photons: A region containing isolated neutral hadron decays, e.g. from hard- 
fragmenting jets, as well as real photons that have a shower-shape fluctuation that fails the tight 
selection. 

• D: nontight, nonisolated photons: A region populated by neutral hadron decays within jets, but 
which have both a small admixture of photons that fail the tight selection and are accompanied by 
a local upward fluctuation of the UE. 

The nontight and nonisolated photons are used to estimate the background from jet events in the signal 
region A. This is appropriate provided there is no correlation between the axes for background photon 
candidates, e.g. that the probability of a neutral hadron decay satisfying the tight or nontight selection cri¬ 
teria is not dependent on whether or not it is isolated. This was studied using a sample of high-pT photon 
candidates from the large sample of Pythia inclusive jet events. Possible correlations, parameterized by 
the Rbkg ratio [21], Rbkg = are taken as a systematic uncertainty, as discussed in 

Sect. 7. 


If there is no leakage of signal from region A to the other nonsignal regions (B, C and D), the double 
sideband approach utilizes the ratio of counts in C to D to extrapolate the measured number of counts in 
region B to correct the measured number of counts in region A, i.e. 


^Sig ^ 


- 


ly^obs 


Nl 


obs 


jw^obs' 


(3) 
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Leakage of signal into the background regions needs to be removed before attempting to extrapolate into 
the signal region. A set of “leakage factors” Ci are calculated to extrapolate the number of signal events 
in region A into the other regions. The leakage factors are calculated using the Pythia simulations in 
intervals of reconstructed photon pj as Ci = where is the number of simulated tight, 

isolated photons. In the 40-80% centrality interval, for Ipl < 1.37 and for 22 < pj < 280 GeV, cb is 
generally less than 0.01, cq ranges from 0.09 to 0.02, and cd is less than 0.003. In the 0-10% centrality 
interval and over the same px range, cb ranges from 0.08 to 0.11, cc ranges from 0.13 to 0.04, and cd is 
0(1%) or less. Except for cb, which reflects the different isolation distributions in peripheral and central 
events, the leakage factors are of similar magnitude to those derived in the pp data analysis [21]. 

Including these factors, and the correlation parameter /?bkg> the formula becomes 




(N^ - ccQ 


(4) 


Equation (4) is solved for the yield of signal photons with /?bkg assumed to be 1.0. The statistical 
uncertainties in the number of signal photons for each centrality, p and px interval are evaluated with 
5 000 pseudo-experiments. Eor each pseudo-experiment, the parameters and are 

sampled from a multinomial distribution with the probabilities given by the observed values divided 
by their sum. The values of Ag®, A^‘®, and A^® used to determine the leakage factors in each 
experiment, are themselves sampled from a Gaussian distribution with the parameters determined by the 
means of the simulated distributions and their statistical uncertainties. Pseudo-experiments where the 
leakage correction is negative are discarded, to exclude trials where the extracted yield is larger than 
A°^^ The standard deviation of the distribution of A^® obtained from the set of pseudo-experiments is 
taken as the statistical uncertainty. 

The purity of the photon sample in the double sideband method is then defined as P = A^‘®/A^’’^ The 
exfracfed values of P are shown in Pig. 4 as a funclion of fransverse momenfum in fhe four measured 
cenfralify infervals and fwo p intervals. In all four cenfralify and bofh p intervals, fhe purify increases 
from abouf 0.5 af fhe lowesf px interval fo 0.9 af fhe highesf px intervals, wifh fypically lower values in 
fhe forward p region. The sfafisfical uncertainly in fhe purify is determined specifically using fhe pseudo- 
experimenfs described above, and by using fhe boundaries defined by fhe highesf and lowesf 16% of fhe 
purify disfribulions lo delermine fhe upper and lower asymmefric error bars. 

Eor kinemafic regions in which fhe number of candidafes in fhe sidebands are small, particularly af fhe 
highesf Px values, fhe populafion of fhose sidebands are re-esfimafed using a dafa-driven approach. Eor 
Ibis, fhe ratio of each sideband (B, C, and D) fo region A as a funclion of px is measured and exfrapolafed 
linearly in 1/px, utilizing all of fhe available dala up lo px = 140 GeV. If should be noted lhal fhe purify 
merely represenls fhe oulcome of fhe sideband sublracfion procedure, and is nol used as an independenl 
correction factor. The several poinls for which Ihis exfrapolafion is ulilized are fhose lo fhe righf of fhe 
verfical doffed line in several of fhe Pig. 4 cenfralify infervals. 

The reconslrucfion efficiency is fhe fraclion of fighl, isolaled phofons malched lo fhe Irulh pholons defined 
above (P^° < 6 GeV), according lo fhe crilerion specified in Seel. 4. The Irue photon px is used in fhe 
numeralor and fhe denominalor, while fhe reconslrucfed p is used in fhe numeralor fo eslimafe fhe very 
small inflow and oufflow of photons in fhe large p infervals used in fhe analysis. The folal efficiency can 
be facforized info fhe producl of Ihree conlribulions: 
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Figure 4: (Color online) Photon purity as a function of collision centrality (left to right) and photon px, for photons 
measured in I 77 I < 1.37 ((a)-(d)) and 1.52 < I;;! < 2.37 ((e)-(h)). The px intervals to the right of the vertical dotted 
line indicated in some bins use the extrapolation method described in the text to account for low event counts in the 
sidebands. 


• Reconstruction efficiency: the probability that a photon is reconstructed with a px greater than 
10 GeV. In the reconstruction algorithm, the losses primarily stem from a subset of photon con¬ 
versions, for which the photon is reconstructed as an electron (“photon to electron leakage”). The 
losses are typically 5% near rj - 0 and increase to about 10% at forward angles, and are found to 
be approximately constant as a function of transverse momentum and centrality. 

• Identification efficiency: the probability that a reconstructed photon passes the tight identification 
selection criteria. 

• Isolation efficiency: the probability that a photon that would be reconstructed and pass the iden¬ 
tification selection criteria, also passes the chosen isolation selection. The large fluctuations from 
the UE in heavy-ion collisions can lead to a photon being found in the nonisolated region. 

Figure 5 shows the total efficiency for each centrality and 77 interval as a function of photon px- The 
primary systematic uncertainties on the efficiency were evaluated by removing the small correction factors 
applied to the simulated shower shapes, and by excluding fragmentation photons from the sample used 
to derive the efficiencies. The contribution from each individual shower-shape selection is small, and so 
the effect on the efficiency is typically small, but the cumulative effect is as large as 10 % in the lowest px 
intervals in the forward p region. Similar correction factors were calculated using the Sherpa simulations, 
and they are found to be consistent with the Pythia calculations in all considered centrality and p regions. 


To account for the residual deviations of the measured photon px from the true px, stemming primarily 
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Figure 5: (Color online) Total photon efficiency as a function of photon px and event centrality averaged over 
I 77 I < 1.37 ((a)-(d)) and 1.52 < |? 7 | < 2.37 ((e)-(h)). Variations of the efficiency from removing the small corrections 
to the simulated shower-shape variable, and from removing fragmentation photons from the simulations are shown 
by dotted and dashed lines, respectively. 


from converted photons treated as unconverted, and from the photon energy resolution, the data are cor¬ 
rected using a bin-by-bin correction technique [21] to generate the correction factors ’li. For each interval 
in centrality and p, a response matrix is formed by correlating the reconstructed px with the truth px for 
truth-matched photons. The projections onto each px interval along the truth axis T{ and the reconstructed 
axis are then constructed for each centrality and p interval and their ratio C\ = T^IRi is formed to calcu¬ 
late the correction in the corresponding px interval. To reduce the effect of statistical fluctuations, the Q 
values were fit to a smooth functional form before applying to the data, with the deviations of the extracted 
correction factors from the fit being generally 0(1 %). In the lowest px interval (22.1 < px < 28 GeV), the 
correction factors deviate from unity by -i-(6-9)% in the central p region and -r(8-13)% in the forward p 
region (the first number for the 40-80% centrality interval, and the second for the 0-10% interval). They 
approach unity rapidly as a function of px and in the highest px interval are -2% in the central p region 
and -1-2% in the forward p region. The reconstructed spectral shapes were compared between simulation 
and data, and were found to agree within statistical uncertainties. Thus, no reweighting of the simulated 
spectrum was performed before calculating the bin-by-bin factors. 

Samples of simulated W and Z bosons decaying to electrons or positrons, based on Powheg [47] interfaced 
to the Pythia8 generator (version 8.175) [48], were used to study the estimated contamination rate relative 
to the total photon rates expected from Jetphox. The raw contamination electron rates were corrected 
using the photon total efficiency. The difference in the extracted cross section of contamination electrons 
between the most peripheral and the most central events was found to be modest. Therefore the centrality 
dependence is neglected and the cross sections calculated for the most central events are used in all 
centrality intervals. Based on this study, it was estimated that the largest background of the W and Z 
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Table 2; Relative systematic uncertainties, expressed as a percentage, on the efficiency-corrected yields for selected 
Pj and centrality intervals in the two 77 intervals. 


n 

|77| < 1.37 

1.52 <|77|< 2.37 

Centrality 

40-80% 

0 - 10 % 

40-80% 

0 - 10 % 

Pt [GeV] 

22-28 

55.6-70 

22-28 

70-88.2 

22-28 

55.6-70 

22-28 

70-88.2 

7 —> e leakage 

1 

1 

1 

1 

1 

1 

1 

1 

Shower shape corn 

3 

2 

5 

3 

6 

2 

9 
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Isolation 

7 

5 

6 

8 

6 

10 

5 

9 

Frag, photons 

< 1 

< 1 

1 

2 

1 

< 1 

2 

2 

Leakage factors 

10 

4 

12 

9 

7 

1 

15 

10 

Nontight criteria 

4 

4 

3 

3 

7 

6 

6 

5 

Rbkg 

21 

7 

13 

6 

20 

4 

15 

11 

Energy scale 

5 

5 

5 

5 

5 

5 

5 

5 

VF/Z contamination 

< 1 

1 

< 1 

1 

1 

1 

1 

1 

Cent, weight 

4 

1 

1 

< 1 

3 

1 

4 

< 1 

77 leakage 

< 1 

< 1 

< 1 

< 1 

2 

1 

2 

2 

Total [%] 

26 

12 

21 

15 

25 

14 

25 

19 


Table 3: Relative systematic uncertainties, expressed as a percentage, on the ratio of the yields in the forward rj 
region and those in the central rj region Rpcr; for selected px and centrality intervals in the two rj intervals. 


Centrality 

1 40-80% 

0 - 10 % 1 

Pt [GeV] 

22-28 

55.6-70 

22-28 

70-88.2 

7 —> e leakage 

1 

1 

1 

1 

Shower shape corn 

3 

0 

4 

0 

Isolation 

9 

9 

4 

4 

Frag, photons 

1 

0 

1 

1 

Nontight criteria 

3 

3 

4 

4 

Leakage factors 

2 

2 

2 

4 

t^bkg 

1 

4 

2 

6 

Energy scale 

7 

7 

7 

7 

VF/Z contamination 

0 

1 

0 

1 

Cent, weight 

1 

0 

4 

1 

77 leakage 

2 

1 

2 

2 

Total [%] 

13 

13 

11 

11 


background is expected in the 35 < px < 44.1 GeV interval with a magnitude of about 8 % in the forward 
pseudorapidity region, and about 5% in the central region. In other bins the correction is smaller, and in 
most bins it is less than 2% in the central 77 region and less than 3% in the forward 77 region. 


7 Systematic uncertainties 


The following systematic uncertainties are accounted for in this analysis. They are broadly classified into 
uncertainties that affect the efficiency, those that affect the yield extraction, and several other additional 
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effects. 


The systematic uncertainties that primarily atfect the total efficiency are: 

• Photon-to-electron leakage: The misidentification of photons as electrons, due to conversions, was 
studied using a sample simulated with extra material, and is found to be less than a 1 % etfect on 
the reconstruction efficiency, since these photons are considered unrecoverable. 

• Shower-shape corrections: To assess the cumulative etfect of the small shower-shape corrections 
applied to mitigate the ditferences between data and simulation, the corrections are removed and 
the ditference in the recalculated yields taken as a conservative systematic uncertainty. This is a 
smaller etfect at higher pj but is as large as 9% at low pT in the forward 77 region. 

• Isolation criteria: To assess the impact of differences between the underlying Ej° distributions in 
data and simulation, several changes in the isolation selection were made. In one case, the cone 
size was changed to A/?iso = 0.4 and the selection enlarged to 10 GeV. In the second, the 
E}^° selection was varied up and down by 2 GeV. Finally, the gap along the E}^° axis between 
regions A/C and B/D was removed. In all of these cases, the selections were similarly adjusted 
in simulation. In general, the variations in the yields show only a weak dependence on px. To 
reduce the effect of statistical fluctuations, the variations as a function fo pj are fit to constants 
over 22 < pt < 44.1 GeV and 44.1 < pj < 140 GeV, and the most signficant variation is applied 
symmetrically to all points in that pj region. If the fit value is consistent with zero, then the 
variation is reduced by half to avoid overcounting the statistical fluctuations. For the foward-central 
ratios, the varations are fit with a single function over 22 < pj < 70 GeV. In several cases, 
changing the isolation selection led to 0 ( 10 %) changes that were clearly consistent with statistical 
fluctuations. In these cases, the variation was reduced to be 5%, similar to the adjacent centrality 
interval. 

The shower leakage corrections were been varied by 1 % of the measured photon px in data, but not 
in simulation, to account for possible defects in the correction. 

• Fragmentation contribution: Excluding the fragmentation photons from the simulation sample has 
typically less than a 2 % effect on the final yields over the full px range. 

The systematic uncertainties that primarily affect the purity of the photon sample in each kinematic and 
centrality interval are: 

• Leakage factors: To test the sensitivity to mismodeling of the shower fluctuations that lead to 
leakage into sideband regions C and D, the leakage factors were conservatively varied up and down 
by 50%. The magnitude is given by the difference between the leakage factors in the 40-80% 
peripheral events, where the underlying event does not cause large extra fluctuations, and the 0 - 
10% most central events. This leads to up to 10% variations at low px, while the effect at higher 
Px is below 5%. 

• Nontight definition: In order to assess the sensitivity to the choice of nontight criteria, which allow 
background into the analysis, the nontight definition was changed from four reversed conditions, to 
five (adding Ws.tot) and two (using just Egide and 0 ) 5 , 3 ). Similar to isolation criteria variations, fits to 
constant values in two px intervals (and one interval for the forward-central ratios) were performed 
to smoothen the bin-to-bin statistical fluctuations. In the central rj interval, the variation is typically 
less than 5% while it is 7% or less in the forward rj interval. 
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• Correlation of tight and isolation axes: The large inclusive-jet Pythia samples were used to study 
possible correlations between the tight selection criteria and the isolation transverse energy. This is 
characterized by calculating 7?bkg for the backgrounds from jets, where the candidate is not matched 
to a truth photon. After integrating over centrality and px, f^bkg was found to vary by about 10% 
in the central rj region and 20% in the forward ij region, albeit with large statistical uncertainties. 
A conservative variation of +20% was propagated through the analysis, which gives up to a 20% 
change at low pj where the purity is lowest, decreasing to typically less than 10% at higher px- 

Uncertainties that pertain to corrections on the energy scale, electron contamination and centrality are 
described here. 

• Energy scale and resolution corrections: The effect of the energy scale and resolution from vari¬ 
ations in material, different energy calibration schemes, and known differences between data and 
simulations in pp collisions are propagated into the bin-by-bin correction factors. The overall vari¬ 
ation from the known sources is typically found to be below 2-3%, and is approximately constant in 
px, but grows at high px in the forward p region. However, the extra-material sample shows a small, 
but systematic, overall shift in the reconstructed energy scale which is approximately independent 
of Px and centrality, but is larger in the forward rj interval. Based on this, an overall uncertainty of 
5% is assigned in the central p region and in the forward region, except in the forward region above 
88.2 GeV, where 7% is assigned. In the ratio, these errors are treated as fully uncorrelated between 
the two p regions. 

• Electron contamination: The contamination from W and Z bosons was estimated to be largest in 
the two Px intervals between 35 and 55.6 GeV, and smaller in the other px intervals. Since the 
calculation does not account for the different expected leakage of the electrons into the different 
sidebands, and since the number of Z bosons in the heavy-ion data is too low to determine this 
fully, 50% of the contamination has been assigned as an uncertainty, leading to a maximum of 4% 
in one px interval in the forward region and smaller in all other intervals. 

• Centrality: The uncertainty on (Taa) for each centrality interval is given in Table 1 and is shared by 
all Px and p intervals for that centrality interval. In addition, the effect of reweighting the simulated 
ECal distribution generally has a less than 2% effect on the final yields, although the impact can 
increase to up to 4% at low px in the forward p interval. 

• p leakage: To address the effect of photons migrating in and out of the large p intervals when 
calculating the efficiency, the true p was also used for the efficiency calculations and was found to 
have a 1-2% overall effect, reaching the larger end of this range in the forward p region. 

Eor the absolute yields, all contributions are added in quadrature. Eor Rpcrj, the systematic variations 
are performed based on the ratio of the forward and central p intervals after each variation, to account 
for correlations between the two p regions. Thus, several of the effects discussed above, particularly the 
influence of the variations in the identification and isolation selection, partially cancel. 

In the central p region, the uncertainties at lower px range from 18% to 26%, and those at higher px range 
from 8% to 16%. In the forward p region, the uncertainties at lower px range from 20% to 26%, and 
those at higher px range from 13% to 19%. Eor the yields, uncertainties for specific centrality, p and px 
ranges are provided in Table 2. Eor the ratio Rfct], the uncertainties at lower px range from 8% to 17% 
and at higher px from 6% to 12%. Uncertainties for specific centrality, p and px ranges are provided in 
Table 2. Eor the ratios, uncertaities for specific centrality and px ranges are provided in Table 3. 
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8 Theoretical predictions 


Jetphox 1.3 is used for NLO pQCD calculations to compare with the fully corrected measurements. 
Jetphox was found to agree well (within 10-15%) with pp from the Tevatron [19, 20] and pp data from 
the LHC [21, 22, 23]. It provides access to a wide range of existing PDF sets and performs calculations 
for direct photon production as well as for photons from fragmentation processes, both using an imple¬ 
mentation of the experimental isolation selection built into the calculations. The primary pp calculations 
shown in this work use the CTEQ6.6 [49] proton PDF, with no nuclear modification, and the BFG 11 
fragmentation functions [50]. They require less than 6 GeV isolation energy in a cone of A/?iso = 0.3 
relative to the photon direction. The effect of hadronization on the final cross secfions was estimated 
using the Pythia 6.423 simulations to be 1% or less, and is neglected in the results shown here. Scale 
uncertainties are estimated by varying the renormalization (//r), factorization (jUp) and fragmentation (/if) 
scales by a factor of two, relative to the baseline result, pr - = pf = Two types of variations 

are performed, a correlated variation of all three scales by a factor of two up and down, as well as an 
independent variation of each scale up and down by a factor of two, leaving the other two scales constant. 
The envelope covered by these variations is typically 12-18%, varying with rj and pj -. PDF uncertainties 
are determined by varying the PDF fit parameters according to 22 eigenvectors in the parameter space, 
and separately keeping track of the upward and downward variations of the final cross secfions. These un- 
cerfainfies are generally less fhan 3% for pj < 100 GeV buf increase fo 6% for pp for px > 140 GeV. The 
impacf of fhe uncerfainfy in fhe sfrong coupling consfanf as{Mz), Aas - ±0.0012, was defermined and 
found fo be small. For fhe yields if varies from ±(l-2)%, decreasing wifh pj . For fhe rafio, if increases 
wifh Pj from 0 fo 2.5%. These errors are nof incorporafed in fhe error bands shown. The calculafions 
were also performed wifh fhe MSTW2008 PDF [51], which yield cross secfions abouf 6% higher for 
|?7l < 1.37 for all calculated pj values. 

To sfudy nuclear effecls, fwo additional calculations are performed. The firsl reweighfs fhe confribufions 
from up and down valence quarks fo accounf for fhe neufrons in fhe colliding lead nuclei, buf wifh no 
affempf af modeling fhe impacf parameter dependence of fhe neufron spatial disfribufions, e.g. due fo a 
neufron skin. This is a reasonable firsf-order approximation for Pb±Pb (bofh wifh A=208) collisions using 
fhe sfandard PDF. The ofher incorporates nuclear modificafions fo fhe nucleon parfon disfribufions using 
fhe EPS09 [1] PDE sef, which are x- and 2^-dependenf modificafions of fhe CTEQ 6.1 PDE, defined 
as ratios of fhe sfandard PDE as a funcfion of x af a hardness scale - 1.69 GeV^ and evolved fo 
fhe relevanf using sfandard DGEAP evolution. The EPS09 modificafions have fheir own sef of 15 
uncertainly eigenvectors, which are used to evaluale 30 variations of fhe cross secfions relalive fo fhe 
defaull sef, which are typically approximately 5%, wifh only a small variafion in pj . 


9 Results 

The per-evenf differenfial phofon yields are calculafed according fo Eq. (2). These are fhen divided by 
(Taa) for comparison wifh fhe Jetphox calculafions. The resulfs are shown as a funcfion of pj in Eig. 6 
and are labulaled in Table 4 for \ r ]\ < 1.37 and in Table 5 for 1.52 < \ r ]\ < 231 . Each panel shows a 
single pseudorapidily inferval, wifh four cenfralily inlervals, each scaled by a factor of 10 relative to each 
ofher. 

The ratios of fhe dafa fo fhe Jetphox pp predicfions are shown in Eig. 7, and are labulaled in Table 6 
for \t]\ < 1.37 and in Table 7 for 1.52 < \q\ < 2.37. In addition, fhe fwo ofher Jetphox calculations 
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Photon [GeV] 



Photon p^ [GeV] 


Figure 6: (Color online) Fully corrected yields of prompt photons in four centrality intervals as a function of px in 
I 77 I < 1.37 (a) and 1.52 < |? 7 | < 2.37 (b) using tight selection, isolation cone size A/?iso = 0.3 and isolation transverse 
energy of less than 6 GeV. Jetphox calculations, for proton-proton collisions and using the same isolation criterion, 
are shown by the yellow bands. Statistical uncertainties are shown by the error bars. Systematic uncertainties on 
the photon yields shown by braces, which are smaller than the markers for some points. The scale uncertainties due 
to (Taa) are tabulated for each bin in Table 1. 


described in the previous section are shown, also divided by the pp results: Pb-i-Pb collisions with no 
nuclear modification (black line), and Pb-i-Pb collisions with EPS09 nuclear modifications (hatched blue 
area). The combined scale and PDF uncertainty on the Jetphox calculations, calculated separately for 
each configuration, are shown as shaded regions. The data are found to agree well with the Jetphox pp 
prediction in all centrality and rj regions, within the stated statistical and systematic uncertainties. They 
are also consistent within uncertainties of the other physics scenarios as well. Thus, the current data are 
not of sufficient precision to address nuclear PDF effects quantitatively. However, it should be noted that 
where the data are more precise, in the central rj region, the PDF modifications implemented in EPS09 
compensate for the suppression at higher px seen in the Pb-i-Pb calculations, giving cross sections similar 
to the pp case. 

The ratios Rpcrj of cross sections between the forward and central rj intervals, are calculated as a function 
of pt for each centrality interval, and are shown in Fig. 8, as well as tabulated in Table 8. Evaluation 
of these ratios leads to the cancellation of several systematic effects on the efficiencies and bin-by-bin 
correction factors, mitigate the effect of the theoretical uncertainties, and fully remove the uncertainty on 
(Taa)- The results are compared to Jetphox calculations for pp (yellow region), Pb-i-Pb (black line with 
grey area) and EPS09 nPDF (grey line with blue area). It is clear that there is some sensitivity to the 
nuclear PDF, primarily through the expected depletion of photon yields in the forward direction expected 
when including the neutron PDF to match the isospin composition of the lead nuclei. While the data are 
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Figure 7: (Color online) Fully corrected normalized yields of prompt photons as a function of px in \n\ < 1-37 
((a)-(d)) and 1.52 < \q\ < 2.37 ((e)-(h)) using tight photon selection, isolation cone size A/?iso = 0.3 and isolation 
transverse energy of less than 6 GeV, divided by Jetphox predictions for pp collisions, which implement the same 
isolation selection. The combined scale and PDF uncertainty on the Jetphox calculation is shown by the grey line 
with yellow area. In addition two other Jetphox calculations are shown, also divided by the pp results: Pb+Pb 
collisions with no nuclear modification (black line with grey area), and Pb+Pb collisions with EPS09 nuclear 
modifications (grey line with blue area). Statistical uncertainties are shown by the bars. Systematic uncertainties 
on the photon yields are combined and shown by the upper and lower braces. The scale uncertainties due only to 
(Taa) are tabulated for each bin in Table 1. 


consistent with all three curves within the statistical and systematic uncertainties, a slight preference for 
the calculations incorporating isospin elfects is observed. 
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Table 6; (7’AA)-scaled prompt photon yields divided by the cross section from pp Jetphox 1.3, for I;;! < 1.37 in four 
centrality intervals as a function of photon p^. 



dN ! dpi! {Txa) ! dcr ! dpi( 5 'e.Temix) \ 

PT [GeV] 

40-80% 

20-40% 

10 - 20 % 

0 - 10 % 

Jetphox 

22-28 

0.95 ± 0.09 ± 0.24 

1.01 ±0.05 ±0.22 

1.15 ±0.04 ±0.20 

1.07 ±0.04 ±0.22 

1 +0.15 
^-0.15 

28-35 

1.04 ± 0.09 ±0.18 

1.08 ±0.05 ±0.17 

1.07 ±0.06 ±0.16 

1.13 ±0.05 ±0.19 

1 +0.14 
^-0.14 

35^4.1 

1.05 ±0.11 ±0.16 

1.08 ±0.06 ±0.14 

1.14 ±0.06 ±0.15 

1.16 ±0.05 ±0.16 

1 +0.13 
^-0.13 

44.1-55.6 

1.10±0.11 ±0.13 

1.06 ±0.07 ±0.12 

1.17 ±0.08 ±0.15 

1.13 ±0.07 ±0.17 

1 +0.15 
^-0.15 

55.6-70 

1.10±0.18±0.13 

1.13±0.10±0.13 

1.05 ±0.10 ±0.12 

1.15±0.11 ±0.18 

1 + 
P P 

70-88.2 

1.33 ±0.23 ±0.14 

1.15±0.18±0.14 

1.23 ±0.13 ±0.15 

1.10±0.13±0.16 

1 + 
P P 

88.2-140 

0.76 ± 0.23 ± 0.08 

1.06 ±0.13 ±0.09 

0.80 ±0.19 ±0.10 

1.05 ±0.13 ±0.12 

1 +0.12 
^-0.12 

140-280 


0.97 ±0.30 ±0.08 

1.22 ±0.43 ±0.12 

1.06 ±0.27 ±0.11 

1+0.15 

^-0.14 


Table 7: (rAA)-scaled prompt photon yields divided by the cross section from pp Jetphox 1.3, for 1.52 < |? 7 | < 2.37 
in four centrality intervals as a function of photon p^. 



dN ! dpi ! {Taa) ! dcr ! dpi (Jetphox) I 

Pi [GeV] 

40-80% 

20-40% 

10 - 20 % 

0 - 10 % 

Jetphox 

22-28 

0.86 ±0.14 ±0.22 

0.90 ±0.07 ±0.19 

0.83 ± 0.08 ± 0.22 

0.91 ±0.07 ±0.23 

1 +0.17 
^-0.17 

28-35 

0.79 ±0.16 ±0.19 

0.89 ±0.09 ±0.18 

0.85 ±0.10 ±0.22 

0.84 ± 0.09 ± 0.22 

1 +0.16 
^-0.16 

35^4.1 

0.74 ±0.20 ±0.17 

0.84 ±0.12 ±0.17 

0.97 ±0.10 ±0.19 

0.78 ±0.14 ±0.16 

1 +0.14 
^-0.14 

44.1-55.6 

0.75 ±0.27 ±0.17 

0.99 ±0.13 ±0.16 

1.16 ±0.08 ±0.16 

1.00 ±0.09 ±0.15 

1 +0.17 
^-0.17 

55.6-70 

0.92 ±0.20 ±0.13 

0.95 ±0.23 ±0.15 

1.24 ±0.13 ±0.17 

0.67 ± 0.24 ± 0.09 

1 + 
P P 

70-88.2 

0.80 ±0.27 ±0.12 

0.76 ±0.20 ±0.10 

1.11 ±0.19±0.14 

0.82 ± 0.20 ±0.16 

1 +0.17 
^-0.17 

88.2-140 

0.64 ± 0.34 ± 0.09 

1.08 ±0.29 ±0.15 

0.97 ±0.26 ±0.12 

0.84 ±0.27 ±0.15 

1 +0.15 
^-0.15 


Table 8 : Results for R^cri, the prompt photon yield in the forward 77 region divided by that in the central rj region, 
as a function of photon pj for four centrality bins ((a)-(d)). Jetphox 1.3 pp calculations are also provided. 



RfCi] 

= dNldpi(1.52 < \ti\ < 2.37)/dN/dpT(\ri\ < 1.37) 


Pi [ GeV] 

40-80% 

20-40% 

10 - 20 % 

0 - 10 % 

Jetphox 

22-28 

0.54 ±0.10 ±0.07 

0.53 ± 0.05 ± 0.05 

0.43 ± 0.05 ± 0.06 

0.52 ± 0.04 ± 0.06 


28-35 

0.45 ±0.10 ±0.07 

0.49 ± 0.05 ± 0.05 

0.47 ± 0.06 ± 0.06 

0.44 ± 0.05 ± 0.06 

- 0.02 

35^4.1 

0.41 ±0.12 ±0.06 

0.45 ± 0.07 ± 0.05 

0.49 ± 0.06 ± 0.06 

0.39 ± 0.07 ± 0.07 

n co+0.03 
*-*•^"- 0.02 

44.1-55.6 

0.38 ±0.14 ±0.06 

0.52 ± 0.08 ± 0.05 

0.55 ± 0.05 ± 0.05 

0.49 ± 0.05 ± 0.05 

'-'•^'^- 0.02 

55.6-70 

0.42 ±0.11 ±0.06 

0.43 ±0.11 ±0.04 

0.60 ± 0.09 ± 0.06 

0.30 ±0.11 ±0.04 

0 51+0.04 

70-88.2 

0.26 ±0.10 ±0.03 

0.29 ± 0.09 ± 0.02 

0.40 ± 0.08 ± 0.03 

0.33 ± 0.09 ± 0.04 

0 44+0 02 
- 0.02 

88.2-140 

0.28 ±0.17 ±0.04 

0.34 ±0.10 ±0.03 

0.40 ±0.14 ±0.04 

0.27 ± 0.09 ± 0.03 

'-'-’-’_ 0.02 
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Photon [GeV] Photon p^ [GeV] Photon p^ [GeV] Photon p^ [GeV] 

Figure 8; (Color online) Fully corrected yields of prompt photons as a function of px in 1-52 < It/I < 2.37 divided by 
that measured in |?/| < 1.37 using the tight photon selection, isolation cone size ARi^o = 0.3 and isolation transverse 
energy of 6 GeV, for four centrality intervals ((a)-(d)). The yield ratio is compared to Jetphox 1.3 predictions 
that implement the same isolation selection, for three different configurations: for pp collisions (grey line with 
yellow area), Pb+Pb collisions with no nuclear modification (black line with grey area), and Pb+Pb collisions with 
EPS09 nuclear modifications (grey line with blue area). Statistical uncertainties are shown by the bars. Systematic 
uncertainties on the photon yields are combined and shown by the braces. 
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10 Conclusions 


In this paper, measured yields of isolated prompt photons in 0.14 nb“' lead-lead collisions recorded 
by the ATLAS detector at the LHC have been presented as a function of collision centrality (in four 
intervals from 40-80% to 0-10%), in two pseudorapidity regions {\t]\ < 1.37 and 1.52 < \r]\ < 2.37) and for 
photon transverse momenta in the range 22 < px < 280 GeV. Photons were reconstructed using the large- 
acceptance, longitudinally segmented electromagnetic calorimeter, after an event-by-event subtraction of 
the average underlying event in each calorimeter layer in small At] intervals. Backgrounds stemming 
from neutral hadrons in jets are suppressed by a tight shower shape selection and by requiring no more 
than 6 GeV transverse energy in a cone of size AR = 0.3 around each photon. The residual hadronic 
background is determined using a double sideband method, and the remaining signal is corrected for 
efficiency and resolution, as well as electron contamination, to arrive at the per-event yield of photons as 
a function of px. in each p and centrality interval. After scaling the yields by the mean nuclear thickness 
(Taa). the px spectrum in each rj and centrality interval is found to agree, within statistical and systematic 
uncertainties, with next-to-leading-order perturbative QCD calculations of proton-proton collisions. The 
data are also compared with calculations that assume the isospin of Pb-i-Pb collisions, as well as the 
calculations for Pb-i-Pb using the EPS09 nuclear modifications of the proton parton distribution functions. 
The ratios of the forward yields to those near midrapidity {Rfcti) are also shown, and are compared to 
the corresponding ratios from Jetphox. The present data are unable to distinguish between the three 
scenarios. However, the overall consistency of the measured yields with Jetphox expectations for all 
centrality intervals demonstrates that photon yields in heavy ion collisions scale as expected with the 
mean nuclear thickness. This provides further support for the interpretation of the clear modification 
of jet yields in Pb-pPb collisions as a function of centrality, relative to those measured in proton-proton 
collisions, as stemming from energy loss in the hot, dense medium [52]. 
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